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ABSTRACT 

We have carried out single and multi-site photometry of the three (3 Cephei stars 
/? and 15 CMa as well as KZ Mus. For the two stars in CMa, we obtained 270 h of 
measurement in the Stromgren uvy and Johnson V filters, while 150 h of time-resolved 
Stromgren uvy photometry was acquired for KZ Mus. All three stars are multi-periodic 
variables, with three {(3 CMa) and four (15 CMa, KZ Mus) independent pulsation 
modes. Two of the mode frequencies of 15 CMa are new discoveries and one of the 
known modes showed amplitude variations over the last 33 years. Taken together, this 
explains the star's diverse behaviour reported in the literature fully. 

Mode identification by means of the amplitude ratios in the different passbands 
suggests one radial mode for each star. In addition, (3 CMa has a dominant £ = 2 mode 
while its third mode is nonradial with unknown £. The nonradial modes of 15 CMa, 
which are £ ^ 3, form an almost equally split triplet that, if physical, would imply that 
we see the star under an inclination angle larger than 55°. The strongest nonradial 
mode of KZ Mus is = 2, followed by the radial mode and a dipole mode. Its weakest 
known mode is nonradial with unknown £, confirming previous mode identifications 
for the star's pulsations. 

The phased light curve for the strongest mode of 15 CMa has a descending branch 
steeper than the rising branch. A stillstand phenomenon during the rise to maximum 
light is indicated. Given the low photometric amplitude of this nonradial mode this is 
at first sight surprising, but it can be explained by the mode's aspect angle. 

Key words: stars: variables: other - stars: early-type - stars: oscillations - stars: in- 
dividual: P CMa - stars: individual: 15 CMa - stars: individual: KZ Mus - techniques: 
photometric 



1 INTRODUCTION 

The P Cephei stars are a group of early B-type stars with 
masses between 9 and 17 Mq that exhibit light, radial ve- 
locity and line-profile variability on time scales between 2-8 
hours (Stankov & Handler 2005). The cause of their vari- 
ability is pulsation in pressure and gravity modes of low 
radial order. About 110 definite members of this group are 
known in our Galaxy (Stankov & Handler 2005, Pigulski 
2005, Handler 2005), but a large number of candidates was 
also revealed in the Large Magellanic Cloud (Kolaczkowski 
et al. 2004). 

One of the major astrophysical applications of the 
(3 Cephei stars is asteroseismology, the determination of the 
interior structure of pulsating stars via theoretical modelling 



of their normal mode spectra. All main-sequence pulsators 
more massive than the Sun that have been studied astero- 
seismically so far, and where unambiguous results on stellar 
structure could be obtained (e.g. see Pamyatnykh, Handler 
& Dziembowski 2004), belong to the /3 Cephei stars. 

The main prerequisite for an asteroseismic study is that 
a number of pulsation modes large enough is known to re- 
move any ambiguities concerning the position of the star in 
the HR diagram, and that these modes are unambiguously 
identified with their pulsational quantum numbers, k, the 
radial overtone of the mode, the spherical degree £, and the 
azimuthal order m. Even among the f3 Cephei stars, such 
objects are not easy to find, which can often be traced to 
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insufficient data that arc incapable of providing the required 
mode frequency resolution and number of mode detections. 

Consequently, it is worthwhile to select potentially 
promising targets and to obtain reasonably largo amounts 
of timc-rcsolvcd photometric and/or spectroscopic data so 
that a reliable assessment of their asteroseismic promise can 
be made. For the present work, we have chosen three stars 
that seemed to justify such an effort. 

The conspicuous naked-eye star (V = 1.98) /3 CMa is 
known as a radial- velocity variable for more than a hundred 
years (see Struve 1950). Shobbrook (1973a) obtained the 
first extensive photometric data set for this star and anal- 
ysed it together with the radial velocities previously pub- 
lished. He found three pulsation periods for the star, and 
suggested that all of them changed somewhat over the years. 
As his measurements were taken in the V filter only, no pho- 
tometric mode identification could be made, but the differ- 
ent light-to-radial-vclocity amplitudes of the two strongest 
modes noticed by Shobbrook (1973a) suggest that they have 
different i. However, we should note that this conclusion is 
based on the assumption that both modes had the same in- 
trinsic amplitude at the different epochs of the photometric 
and spectroscopic measurements. 

15 CMa {V — 4.80) shows notoriously difficult variabil- 
ity according to the literature. This began with the first 
claim of its radial velocity variability (Campbell 1911) that 
was later retracted (Moore 1936). Among the successive 
studios of its variability, the photometric investigation by 
Shobbrook (1973b) was the most extensive. He identified 
the period of the dominant variation in his light curves, and 
noted its non-sinusoidal shape, but could not push the anal- 
ysis further. Heynderickx (1992) suggested the presence of 
at least two pulsation modes, which wore later identified 
with i = 2 and 0, respectively, on the basis of multicolour 
photometry (Heynderickx, Waelkens & Smeyers 1994). 

KZ Mus {V = 9.1) has been studied during a two-site 
photometric effort by Handler et al. (2003). These authors 
also review the (short) history of the investigations of its 
variability. They found four independent modes in the stel- 
lar light variations, where the strongest of these have i = 2, 
and 1, respectively, and the radial mode is either the fun- 
damental or the first overtone. The presence of additional 
pulsation modes was suspected by Handler et al. (2003), 
who also pointed out the importance of an identification of 
the spherical degree of the fourth known mode. 



2 OBSERVATIONS AND REDUCTIONS 

We acquired time-resolved differential photoelectric pho- 
tometry of all the three targets. The bright stars (3 and 
15 CMa were measured with respect to the comparison stars 
HR 2271 (V = 5.8, B3H/HI) and 17 CMa (V = 5.8, A2V), 
that were observed alternating with the targets. No evi- 
dence for variability of the latter two stars within a limit 
of 1.5 mmag was found. 

This group of four stars was observed from three ob- 
servatories in the time span of November 2004 to Febru- 
ary 2005. Thirty-five nights of Strorngrcn uvy data were 
obtained at Siding Spring Observatory, Australia, with the 
0.6-m telescope. Eighteen nights were acquired with the 0.5- 
m telescope at the South African Astronomical Observatory. 



The Johnson V filter was used for the first five nights; the re- 
maining 13 nights consist of Stromgren uvy data, but the u 
filter was not used for /? CMa at this site. Finally, ten nights 
of measurement were performed with the 0.75-m Automated 
Photometric Telescope T6 at Fairborn Observatory, USA, 
again using Stromgren uvy filters. Neutral density filters 
were used in all cases to avoid photomultiplier destruction, 
owing to the extreme brightness of /3 CMa; the same neu- 
tral density filters were applied to all stars at each site. The 
total of 63 nights of observation has a time base of 111.95 d, 
corresponding to 272.7 h of measurements of /3 CMa and 
265.4 h of measurements for 15 CMa, respectively. 

KZ Mus was measured with the same observational 
technique. The comparison stars wore HD 110022 (V = 7.8, 
B8 HI/IV), which was replaced by HD 109082 {V = 8.1, 
B9IV) after the first few nights when we found evidence for 
variability of the former object. The second comparison star 
was HD 111876 {V = 7.1, B9.5/A0V). HD 109082 was con- 
stant during the observations. However, HD 111876 showed 
a small (~ 0.01 mag iu all three filters) mean magnitude 
drift during the total time span of observations. 

This group of stars was observed from March to June 
2005, but with the 0.6-m telescope at Siding Spring Obser- 
vatory only. Data were acquired during 27 nights, giving a 
total time span of 98.9 d, and 154.3 h of measurement re- 
sulted. Again the Stromgren uvy filters were used, but no 
neutral density filters had to be employed. 

Data reduction comprised a correction for coincidence 
losses, subtraction of sky background and extinction deter- 
mination and correction. Whenever possible, nightly extinc- 
tion coefficients were determined with the classical Bouguer 
method (fitting a straight line to a plot of magnitude vs. air 
mass) from the measurements of the two comparison stars. 
For the remaining nights, mean coefficients from the adjoin- 
ing nights were used. 

The magnitudes of the comparison stars were adjusted 
to their mean differences, and were then combined into a 
curve that was assumed to reflect the effects of sky trans- 
parency and detector sensitivity changes only. Subsequently, 
these combined time series were binned into intervals that 
would allow good compensation for the above-mentioned 
non-intrinsic variations in the target star time series and 
were subtracted from the measurements of the variables. 
The binning minimises the introduction of noise in the dif- 
ferential light curves of the targets. 

The timings for the differential light curves were lie- 
liocentrically corrected as the next step. Finally, the pho- 
tometric zeropoints of the different instruments were com- 
pared between the different sites and were adjusted if nec- 
essary. Measurements in the Stromgren y and Johnson V 
filters were treated as equivalent due to the similar effective 
wavelengths of these filters, and were analysed jointly af- 
ter a correction for differential colour extinction in V. Some 
example light curves of our three targets together with mul- 
tifrequency solutions to be discussed in the following are 
shown in Figs. 1, 2 and 3, respectively. 

3 FREQUENCY ANALYSIS 

Our frequency analysis was mainly performed with the pro- 
gram Period98 (Sperl 1998). This package applies single- 
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Figure 1. Some of our observed light curves of /3 CMa. Plus 
signs are data in the Stromgren u filter, opem circles arc oiu- 
V measurements and filled circles represent Stromgren y data. 
The full line is a fit composed of all the periodicities detected 
in the light curves (Table 2). The amount of data displayed is 
approximately one third of the total. 
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Figure 2. Some of our observed light curves of 15 CMa. Plus 
signs are data in the Stromgren u filter, open circles are our v 
measurements and filled circles represent Stromgren y data. The 
full line is a fit composed of all the periodicities detected in the 
light curves (Table 3). The amount of data displayed is approxi- 
mately one third of the total. 



frequency power spectrum analysis and simultaneous multi- 
frequency sine-wave fitting. It also includes advanced op- 
tions such as the calculation of optimal light-curve fits for 
multi-periodic signals including harmonic, combination, and 
equally spaced frequencies. Our analysis will require some of 
these features. 

3.1 P CMa 

We started by computing the Fourier spectral window of the 
u filter data. It was calculated as the Fourier transform of a 
single noise- free sinusoid with a frequency of 3.979 cd~^ (the 
strongest pulsational signal of /? CMa) and an amplitude of 
12 mmag, sampled in the same way as our measurements. 
The upper panel of Fig. 4 contains the result. The alias 
structures in this window are sufficiently suppressed to allow 
easy and unique determinations of the frequencies of the 
stellar light variations. 

Wc proceeded to compute the amplitude spectra of the 
data themselves (second panel of Fig. 4). The signal desig- 
nated /i dominates. We prewhitened it by subtracting a syn- 
thetic sinusoidal light curve with a frequency, amplitude and 



phase that yielded the smallest possible residual variance. 
We then computed the amplitude spectrum of the residual 
light curve, which is shown in the third panel of Fig. 4. 

A second signal (/2) can clearly be seen in this graph. 
Prewhitening it leaves no further significant periodicities in 
the u data alone. At this stage, we prewhitened the data 
from the v and y filters by the /i and /2 frequencies, then 
combined these residuals with those from the u filter data 
divided by 1.2. This 1.2 factor scales the u data to the same 
amplitude as that found for /i in the v and y data. This 
procedure is valid because the light curves of (3 Cephei stars 
have the same phase at all optical wavelengths. 

The residual amplitude spectrum of the uvy data com- 
bined in this way is shown in the lowest panel of Fig. 4. We 
note the presence of another signal, fs, and include it in our 
frequency solution for /3 CMa. Prewhitening by this varia- 
tion as well, we find no further significant periodicity in our 
combined residual amplitude spectrum. 

Wc consider an independent peak statistically signifi- 
cant if it exceeds an amplitude signal-to-noise ratio of 4 in 
the periodogram; combination signals must satisfy S/N > 
3.5 to be regarded as significant (see Breger et al. 1993, 
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Figure 3. Some of our observed light curves of KZ Mus. Plus 
signs arc data in the Stromgren u filter, opcm circles arc our 
V measurements and filled circles represent Stromgren y data. 
The full line is a fit composed of all the periodicities detected 
in the light curves (Table 4). The amount of data displayed is 
approximately half the total. 



1999 for a more in-depth discussion of this criterion). The 
noise level was calculated as the average amplitude in a 
5 cd"'^ interval centred on the frequency of interest. 

We are now in a position to determine a final multifre- 
qucncy solution for our light curves of B CMa. To this end, 
we first computed a weighted average for the frequencies 
derived from the individual filters' data, where the weight 
corresponded to the S/N of the signals. The results were 
adopted as our final frequency values. We then fitted these 
frequencies to the uvy data and determined their amplitudes 
and phases (that were, as assumed previously, found to be 
the same for all filters within the errors). The final solution 
is listed in the upper part of Table 1. 

The error estimates we give are the formal values for 
the amplitudes (Montgomery & O'Douoghuc 1999). We note 
that such formal error bars axe expected to underestimate 
the real errors by about a factor of 2 (Handler et al. 2000, 
Jerzykiewicz et al. 2005). For the frequencies themselves, 
we computed both the formal errors, and the deviation of 
the mean of the values from the individual filters; we then 
adopted the larger of the two values. For frequencies /i and 
/a the two results agreed quite well (within 10 per cent). 
However, the error of the moan was a factor of 2 larger for /2. 

This result can be easily understood: within the time 
span of our data, /2 is not fully resolved from a possible 
variation of 4 cycles per sidereal day. Such a variation could 
be present in our light curves due to residual (colour) ex- 
tinction effects, despite our careful data reduction. This in- 
terpretation is supported by the observation that the value 
of /2 is closest to 4 cycles per sidereal day in the u band 
where possible residual extinction effects are expected to be 
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Figure 4. Amplitude spectra of /3 CMa. The uppermost panel 
shows the spectral window of the data, followed by the peri- 
odogram of the data. Successive prewhitening steps are shown 
in the following panels; note their different ordinate scales. 



Table 1. Multifrequcncy solutions for photometric data of (3 
CMa. Error estimates (following Montgomery & O'Donoghue 
1999) on the amplitudes are ± 0.4 mmag in u and ± 0.3 mmag 
in V and y, respectively. The S/N ratio was computed following 
Breger et al. (1993) and is quoted for the v filter. 



ID 


Freq. 


u Amp. V Amp. 


y Amp. 


S/N 




(cd-i ) 


(mmag) (mmag) 


(mmag) 




New photometry 


/l 


3.9793 ± 0.0001 


12.0 10.6 


10.0 


25.5 


/2 


3.9994 ± 0.0007 


3.8 1.7 


1.3 


4.1 


h 


4.1857 ± 0.0007 


1.7 1.7 


1.4 


4.2 


Archival data 


/i 


3.9792 ± 0.0001 




10.3 (V) 


28.5 


/2 


3.9999 ± 0.0007 




1.9 (V) 


5.2 


/3 


4.1821 ± 0.0009 




1.4 (V) 


40 
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largest. We must therefore be careful in the interpretation 
of the frequency and uvy amplitudes of /2. 



3.1.1 Reanalysis of archival data 

One of us (RRS) has digitised his archival V filter photome- 
try of l3 and 15 CMa taken in the early 1970s. We can there- 
fore analyse these measurements with our frequency analysis 
methods, allowing a direct comparison of the results. 

For f3 CMa we confirm the findings by Shobbrook 
(1973a, sec lower part of Table 1): wc obtain the same three 
frequencies as he did. Given the probable underestimation 
of the real errors of the frequencies and amplitudes (espe- 
cially for mode /2), we are reluctant to claim the presence 
of amplitude and/or frequency variability for the star. 

On the other hand, wc note that there sccins to bo some 
low-frequency variability, for which we are unable to deter- 
mine a period, remaining in the residuals, similar to what 
can be discerned in the lowest panel of Fig. 4. We find this 
behaviour in all independent data sets wc have for /3 CMa, 
but not in the differential light curve of the two comparison 
stars. Because we can exclude an extinction effect for its ori- 
gin, we are inclined to believe that this variation is intrinsic 
to our target star. 

We attempted to obtain an improved frequency solu- 
tion and to check for amplitude and frequency variability of 
the stellar pulsation modes by analysing the combined vi- 
sual photometric data by Shobbrook (1973a), Balona et al. 
(1996) and the HIPPARCOS satellite (ESA 1997) together 
with our measurements. Regrettably, the amount and tem- 
poral separation of the data render this investigation incon- 
clusive: we derive a best-fitting frequency for the strongest 
mode of 3.9793334 ± 0.0000009 cd"^ that is slightly variable 
in amplitude and phase, but wc cannot be certain whether 
this is a numerical artifact due to beating with the close 
mode at 3.9994 cd"'^ when unresolved in certain subsets of 
data. 



3.2 15 CMa 

The frequency analysis for this star was carried out in a 
similar way as that for /? CMa and is illustrated in Fig. 5. 
This time we use the v data for display as they resulted 
in the light curves of best S/N. In the case of 15 CMa, 
the spectral window had a single noise-free sinusoid with a 
frequency of 5.419 cd"'^ and an amplitude of 5.4 mmag as 
input (upper panel of Fig. 5). The amplitude spectrum of 
the data themselves (second panel) already makes it clear 
that more than one frequency is present in the light curves. 
Prewhitening the strongest signal (/i) leaves two conspicu- 
ous peaJis in the residual amplitude spectrum. One of these 
is the first harmonic of /i . Further prewhitening with these 
three signals indicates the presence of another periodicity, 
and prewhitening that one as well reveals a fourth inde- 
pendent frequency in the light curves. Pushing the analysis 
further, also by combining the residuals for the u, v and y 
data, as we did for f3 CMa, did not result in the detection 
of additional periodicities. 

The final multifrequency solution for 15 CMa was deter- 
mined in the same way as for (3 CMa: by computing S/N- 
weighted averages of the frequencies determined from the 
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Figure 5. Amplitude spectra of 15 CMa. The uppermost panel 
shows the spectral window of the data, followed by the peri- 
odogram of the data. Successive prewhitening steps are shown 
in the following panels; note their different ordinate scales. 



individual filters' data. Wc note that, for reasons of better 
accuracy, we fixed the frequency of the harmonic signal to 
exactly twice the value of the parent mode with Period 98. 

Wc then recomputed the uvy amplitudes with these 
"best" frequencies and list the result in Table 2. Again, the 
adopted error estimates on the frequencies correspond to the 
larger of the two values derived from analytic formulae and 
rms deviations of the mean of the frequencies. 



3.2.1 Reanalysis of archival data 

We analysed the V measurements of 15 CMa by Shobbrook 
(1973b, also digitised by RRS), who also detected our signal 
/i and its first harmonic. However, he could not continue to 
search for further frequencies because the residual amplitude 
spectrum was too complicated. 

Thanks to our multi-site measurements, wc know four 
pulsation frequencies of the star without aliasing ambigui- 
ties, and we can use this information in our reanalysis. This 
helps us to solve the puzzle that the data by Shobbrook 
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Table 2. Multifrequency solution for time-resolved photomet- 
ric data of 15 CMa. Error estimates (following Montgomery & 
O'Donoghue 1999) on the amplitudes are ± 0.3 mmag in u and ± 
0.2 mmag in v and y, respectively. The S/N ratio was computed 
following Breger et al. (1993, 1999) and is quoted for the v data. 
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ID 


Freq. 


u Amp. 


V Amp. 


y Amp. 


S/N 




(cd-i ) 


(mmag) 


(mmag) 


(mmag) 




New photometry 


fl 


5.4187 ± 0.0002 


6.2 


5.4 


4.6 


16.0 


/2 


5.1831 ± 0.0003 


5.4 


3.9 


2.9 


11.2 


2/l 


10.8374 


2.4 


2.1 


2.1 


7.2 


/3 


5.5212 ± 0.0008 


2.0 


1.6 


1.3 


4.7 


/4 


5.3085 ± 0.0012 


1.6 


1.4 


1.1 


4.2 



(1973b) provided: we first detect /i and its liarmonic. Tlie 
three other frequencies detected in onr new measurements 
are also present in the older photometry, and all three are 
required to explain the observations, as proven by different 
prewhitening trials. We list the result of this analysis in the 
lower part of Table 2. We point out that /2 has a S/N below 
4, but since it is known to be an intrinsic signal from our 
new measurements, we only require S/N > 3.5 to accept 
this frequency as real in the old data. 

In this case, /2 can be safely claimed as having varied in 
amplitude, whereas no such statement can be made for the 
other signals. In this context it is interesting to revisit the 
discussion of the Geneva data by Heynderickx (1992), who 
also observed 15 CMa. He detected our /i and its harmonic, 
but then encountered two closely spaced signals around /2, 
where the frequency difference of these signals is close to 
the inverse time span of his data set. Consequently, it can 
be suspected that these two apparent signals near /2 could 
in fact be due to a single mode with amplitude variability 
during Heynderickx' observations. Wc note that an attempt 
to refine the pulsation frequencies of the star by analysing 
the measurements by Shobbrook (1973b) and the HIPPAR- 
COS satellite (ESA 1997) together with ours did not bear 
fruit due to aliasing problems. 

We also note that, just as for /3 CMa, some slow aperi- 
odic variability seems to be left in the residuals of both the 
archival and the new measurements. This slow variability is 
not correlated with the one of /3 CMa (neither in the old nor 
in the new data), which is further support for the idea that 
it is intrinsic to the stars and is not an instrumental effect. 



3.2.2 The light curve shape of fi 

The amplitude of the harmonic frequency of f i is unusually 
high given the amplitude of the independent mode to which 
it is connected. We have therefore examined the shape of 
the light curve due to this mode. We fitted the five known 
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Figure 6. Phase diagram of light curves of 15 CMa with respect 
to fl (full dots), compared to a fit computed with /i and 2/i. 
Two cycles of the variation are shown. The rising branch of this 
light variation is flatter as the descending branch, and there seems 
evidence for a stillstand before light maximum. 



frequencies of 15 CMa to the archival data and to our new 
measurements and then prewhitened the fit composed of the 
resulting parameters for /2, /a and /i. We then phased the 
residuals with /i and show the results in Fig. 6. 

The phase diagrams obtained in this way appear odd 
for a pulsating star: the rising branch of the light curve is 
less steep than the descending branch. When computing a fit 
to these phase diagrams composed of /i and 2/i, or when 
averaging the two phase diagrams, it even seems possible 
that there is a "stillstand" phenomenon in the rising branch. 
We would however like to see more data on the star before 
coming to a definite conclusion on this matter. 

It is interesting to note that the scatter in the two phase 
diagrams axe very similar: although smaller in number, the 
precision of the measurements taken in the 1970s is slightly 
better than that of our present data. 



3.3 KZ Mus 

The frequency analysis for this star was carried out in a 
similar fashion as for the others. We determined individual 
frequencies step-by-step and found four independent modes 
(see Fig. 7) plus three combination frequencies. Again, the 
frequencies were first determined for each individual filter, 
using Period 98 to fix the frequencies of combination sig- 
nals to the exact sums predicted by those of parent modes. 
Then, S/iV-weighted averages of the frequencies determined 
from the individual filters' data were computed and adopted 
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Figure 7. Amplitude spectra of KZ Mus. Tlie uppermost panel 
shows the spectral window of the data, followed by the peri- 
odogram of the data. Successive prewhitening steps are shown 
in the following panels; note their different ordinate scales. 



Table 3. Multifrequency solution for our new time-resolved pho- 
tometry of KZ Mus. Error estimates (following Montgomery & 
O'Donoghue 1999) on the amplitudes are ± 0.4 mmag in u and 
V, and ± 0.3 mmag in y. The S/N ratio was computed following 
Breger et al. (1993, 1999) and is quoted for the u filter data. 



ID 


Freq. 
(cd-1) 


u Amp. 
(mmag) 


V Amp. 
(mmag) 


y Amp. 
(mmag) 


S/N 


/l 


6.86400 ± 0.00004 


47.2 


41.0 


38.4 


85.4 


f2 


5.95059 ± 0.00007 


35.0 


20.6 


17.3 


63.4 


S3 


6.18774 ± 0.00013 


16.1 


12.4 


11.6 


29.3 


/4 


5.7094 ± 0.0005 


4.0 


3.2 


3.6 


7.1 


/l + /2 


11.81459 


1.7 


1.7 


1.7 


3.2 


/l + /3 


12.05174 


2.3 


1.8 


1.7 


4.3 


2/2 


11.90118 


2.2 


1.6 


1.8 


4.1 



present in the older measurements, but we were not con- 
vinced about its reality during the previous analysis. 

On the other hand, the combination signal /2 — /i de- 
tected by Handler et al. (2003) is not present in our new mea- 
surements. This is likely due to the increased low-frequency 
noise in the new data caused by the variability of the com- 
parison star HD 111876, which we were not able to remove 
satisfactorily. 



3.3.1 Refining the determination of the pulsation 
frequencies 

Although the photometry by Handler et al. (2003) was taken 
three years before the measurements presented in this pa- 
per, the frequency determinations of both works are accurate 
enough that we can link the previous V and the present y 
data to obtain more accurate pulsation frequencies. The ef- 
fective wavelengths of these two bands are also sufficiently 
similar to render such an approach valid. 

Consequently, we performed a joint frequency analysis 
with Period 98, where we could also reveal the presence 
of one more combination signal. The result is given in Tar 
ble 4. We have confirmed that the frequencies of the three 
strongest modes and their combination frequencies are alias- 
free by including the HIPPARCOS photometry of the star 
(ESA 1997) as well. We could not make a similar check for 
mode fi due to its low amplitude. In any case, the resid- 
ual amplitude spectrum after prewhitening our joint solu- 
tion contains no evidence that measurable amplitude or fre- 
quency variations have occurred within the 3.4-y time haae 
of our observations. We have also determined the Stromgren 
V amplitudes given this frequency solution for the combined 
2002 and 2005 data. However, we could not do the same for 
the Johnson U and Stromgren u bands because their wave- 
length responses are sufficiently different that the calculated 
amplitudes of the strongest modes do not correspond within 
the errors. 



as our final frequency values. Our multifrequency solution 
obtained this way is listed in Table 3. 

Comparing our results to those by Handler ot al. (2003) , 
we cannot find significant evidence for amplitude or fre- 
quency variability, but we do note the occurrence of a 
new combination frequency, /i -h /s- This signal is actually 



4 MODE IDENTIFICATION 

The spherical degree (. of the pulsation modes of our three 
target stars can be identified by comparing their observed 
uvy amplitudes with those predicted by theoretical models. 
Required input for these models are the mass and effective 
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Table 4. Multifrequency solution for the combined ground-based 
time-resolved V and y, as well as v photometry of KZ Mus. The 
error estimates (Montgomery & O'Donoghue 1999) on the ampli- 
tudes are ±0.2 mmag. The S/N ratio was computed following 
Breger et al. (1993, 1999). 



ID 


Frequency 


V Amp. 


V OT y Amp. 


S/N 




(cd-l) 


(mmag) 


(mmag) 




fl 


6.864016 ± 0.000006 


41.2 


38.4 


109.5 


f2 


6.950693 ± 0.000011 


20.4 


16.5 


47.2 


/3 


6.187472 ± 0.000016 


12.3 


10.9 


31.4 


/4 


5.70935 ± 0.00006 


3.1 


3.0 


8.4 


fl + /2 


11.814709 


1.9 


1.9 


6.2 


2/2 


11.901386 


1.5 


1.5 


5.0 


fl + fs 


12.061489 


1.3 


1.3 


4.4 


f2 + fs 


12.138166 


1.2 


1.2 


4.0 



temperature of the stars, from which their luminosity obvi- 
ously follows as well. 

For KZ Mus, we simply adopt the values given by Han- 
dler et al. (2003), who derived T^s = 26000 ± 700 K and 
log L/Lq = 4.22 ± 0.20. For the two stars in Canis Major, 
the calibrations by Crawford (1978) applied to the stan- 
dard Stromgren photometry by Crawford, Barnes & Gol- 
son (1970) results in E{b - y) = 0.027 for /3 CMa and 
£^(6 ~y) = 0.031 for 15 CMa, respectively. The Stromgren 
system calibration by Napiwotzki, Schonberner & Wenske 
(1993) then yields T^s = 25600 ± 1000 K for (3 CMa and 
Teff = 26300 ± 1000 K for 15 CMa. 

With the Geneva colour indices of the two stars ob- 
tained from the Lausanne Photometric data base, the cal- 
ibrations by Kiinzh et al. (1997) provide Tcff = 25100 ± 
1100 K and log g = 3.3 ± 0.6 for /3 CMa, and Tes = 
26000 ± 1500 K and log 5 = 3.8 ± 0.6 for 15 CMa. 

The analysis of lUE spectra of a number of /? Cephei 
stars by Niemczura & Daszynska-Daszkicwicz (2005) results 
in Teff = 24700 ± 800 K and log 5=3.74 for 13 CMa and 
Teff = 24600 ± 1000 K and log g=3.81 for 15 CMa. 

The HIPPARCOS parallax of CMa is 6.53 ±0.66 mas. 
With 1^=1.98 and the reddening as determined before, we 
obtain M^, = -4.1±0.2. The parallax of 15 CMa (2.02±0.70 
mas) is too inaccurate to be useful. However, 15 CMa is a 
member of the young open cluster Collinder 121, whose dis- 
tance was determined as 592 ±28 pc (de Zeeuw et al. 1999), 
which is in reasonable agreement with the usually quoted 
distance to the cluster in the pre-HIPPARCOS era, 630 pc 
(Feinstein 1967). Using the most recent distance to Collinder 
121, we determine My = —4.2 ± 0.2 for this star, given that 
V^=4.80 and again adopting the amount of reddening sug- 
gested by Stromgren photometry. 

Summarising the temperature determinations quoted 
before, we arrive at final values of T^g = 25100 ± 1200 K 
for (3 CMa and 25600 ± 1200 K for 15 CMa. According 
to Flower (1996), these effective temperatures correspond 
to bolometric corrections of BC = —2.2 ± 0.3 mag and 
BC = —2.3 ±0.4 mag, respectively. Consequently, the bolo- 
metric luminosities of the two stars are —6.3 ± 0.4 mag {/3 
CMa) and -6.5 ± 0.4 mag (15 CMa), respectively. 

Figure 8 shows the positions of our three targets in a 
theoretical HR diagram. The two stars in Canis Major have 
very similar properties there, indicating that both are mas- 
sive and are approaching the end of their main sequence life. 




Figure 8. The positions of our three targets in the theoretical HR 
diagram. The filled circle denotes KZ Mus, the star symbol is for 
fd CMa, and the open circle represents 15 CMa. Error estimates 
on the effective temperatures and luminosities of the stars are 
indicated. Some stellar evolutionary tracks, for a metal abundance 
of Z = 0.02, labelled with their masses (full lines) are included for 
comparison. We also show the theoretical borders of the f3 Cephei 
instability strip (Pamyatnykh 1999, dashed-dotted line). 

This is somewhat in contrast with their distinctly different 
pulsation periods. KZ Mus appears to just have entered the 
instability strip and has lower mass than the two other pul- 

sators. 

For the purpose of mode identification we computed 
theoretical photometric uvy amplitudes for the ^ ^ ^ 4 
modes of models with masses between 13.0 and 16.0 Mq 
in steps of I.OMq, and in a temperature range of 4.38 ^ 
log Toff ^ 4.43 for f3 and 15 CMa. For KZ Mus, we investi- 
gated a range of 11.5 to 14.0 Mq in steps of 0.5 Mq, and 
4.395 < log Teff ^ 4.435, therefore applying a generous in- 
terpretation of the derived error bars. All models had an 
assumed metallicity oi Z = 0.02. This approach is similar to 
that by Balona & Evers (1999). 

Theoretical mode frequencies between 3.7 and 4.5 
cd~^ were considered for f3 CMa, between 4.8 and 5.9 
cd"^ for 15 CMa and between 5.4 and 6.5 cd"^ for KZ Mus. 
The frequency ranges chosen are somewhat larger than the 
ones excited in the stars to allow for some rotational m- 
mode splitting. We compare the theoretical photometric am- 
plitude ratios to the observed ones in Figs. 9 to 11 for all 
independent modes of all three stars. 

Starting with /3 CMa, it is clear that its strongest mode 
has a spherical degree of ^ = 2. The second mode seems ra- 
dial, although its measured amplitude vs. wavelength depen- 
dence is considerably steeper than theoretically predicted. 
However, we believe that this is caused by problems with 
differential colour extinction, which artificially increases the 
measured u amplitude, as discussed in Sect. 3. In any case, 
no mode identification other than £ = is possible for this 



Photometric studies of (3 CMa, 15 CMa and KZ Mus 9 



0.4 



-1 — I — I — r 



] — I — I — I — I — I — I — I — I — I — I — r- 




_L 



1.2 



-p-i— I— I— I— r^r-i— I— I— p-i— I— I— I— p 

/S CMa. 3.9994 c/d 




3500 4000 4500 5000 
Wavelength (A) 



5500 



Figure 9. Mode identifications for /3 CMa from a comparison 
of observed and theoretical uvy amplitude ratios, normalised to 
unity at u. The filled circles with error bars axe the observed 
amplitude ratios. The full lines are theoretical predictions for ra- 
dial modes, the dashed lines for dipole modes, the dashed-dotted 
lines for quadrupole modes, the dotted lines for octupole = 3) 
modes, and the dashcd-dot-dot-dottcd lines are for £ = 4. The 
thin error bars denote the uncertainties in the theoretical ampli- 
tude ratios. 



mode. The lowest-amplitude pulsation mode is nonradial, 
but its spherical degree cannot be constrained due to its 
small amplitude. 

Continuing with 15 CMa, we see that the strongest 
mode is nonradial with £ = 1, 2 or 3, but no distinction be- 
tween these hypotheses can be made. The second strongest 
mode is most likely radial, and no mode identification is 
possible for the two weakest modes. However, due to their 
proximity to the radial mode in frequency, they must be 
nonradial. 

Turning to KZ Mus, one immediately sees that its 
strongest mode \s I = 2. The second strongest mode is sug- 
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Figure 10. Mode identifications for 15 CMa from a comparison 
of observed and theoretical uvy amplitude ratios, normalised to 
unity at u. The filled circles with error bars are the observed am- 
plitude ratios. The full lines are theoretical predictions for radial 
modes, the dashed lines for dipole modes, the dashed-dotted lines 
for quadrupole modes, the dotted lines for f = 3 modes, and the 
dashed-dot-dot-dotted lines are for I = A. The thin error bars 
denote the uncertainties in the theoretical amplitude ratios. 
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Figure 11. Mode identifications for KZ Mus from a comparison 
of observed and theoretical uvy amplitude ratios, normalised to 
unity at u. The filled circles with error bars are the observed am- 
plitude ratios. The full lines arc theoretical predictions for radial 
modes, the dashed lines for dipolc modes, the dashed-dotted lines 
for quadrupole modes, the dotted lines for octuijole modes, and 
the dashed-dot-dot-dotted lines are for £ = 4. The thin error bars 
denote the uncertainties in the theoretical amplitude ratios. 



gested to be radial, although again the observed amplitude 
ratios indicate a steeper wavelength dependence of the am- 
plitudes than theoretically predicted. The third mode of KZ 
Mus is a dipole, £ = 1, and the fourth is nonradial. These re- 
sults arc in perfect agreement with the mode identifications 
by Handler et al. (2003). 



5 DISCUSSION AND CONCLUSIONS 

We have observed the three /3 Cephei stars because we aimed 
at new and improved mode identifications and because we 
wanted to check whether these objects would justify a larger 
observational effort for asteroseismic purposes. We have rea- 
sonably succeeded in providing the mode identifications. 

One point that appears problematic, however, is that 
the theoretically predicted u/v and u/y amplitude ratios for 
radial modes are smaller than the measured ones for l3 CMa 
and KZ Mus. Whereas the reason for this discrepancy may 
be incorrect differential colour extinction coefficients for /3 
CMa that directly affect the measured amplitudes of the 
3.9994cd"^ mode, we cannot invoke this explanation for the 
S.QSlcd^'^ mode of KZ Mus as its frequency is too differ- 
ent from an integer multiple of one cycle per sidereal day. 
In addition, the theoretical and observed UvBV amplitudes 
derived by Handler et al. (2003) arc in perfect agreement for 
a radial mode in this star. We note that we have used the 
same differential colour extinction correction for l3 and 15 
CMa, which was derived from the 15 CMa data only. These 
two stars were measured simultaneously, and have very sim- 
ilar colours, which is why this procedure was soundest. 

Returning to the mode identifications of the individual 
stars, the two strongest modes of /3 CMa are a quadrupole 
and likely a radial mode. Given the effective temperature, 
luminosity and thus mass inferred in Sect. 4, we can com- 
pute the pulsation "constant" for such a radial mode. The 
result of this calculation is Q = 0.035 ± 0.009 d, indicating 
that this mode may be the radial fundamental, although the 
possibility that it is the first overtone cannot be ruled out 
this way. 

We were also able to locate a radial mode in the pulsa- 
tion spectrum of 15 CMa, whose pulsation constant would 
be Q = 0.027 ± 0.007 d. This suggests that we deal with the 
first radial overtone in this case, although the fundamen- 
tal mode or the second overtone can also not be excluded. 
We note, however, that no convincing detection of a radial 
mode in a /3 Cephei star other than the fundamental has 
been made to date. 

Three of the nonradial pulsation modes of 15 CMa are 
almost equally spaced in frequency. Regrettably, we cannot 
be sure whether these modes are £ = 1, 2 or 3. Independent 
of the £ value, we can still assume that these modes do cor- 
respond to rotationally split components of the same mode 
to check whether this leads to reasonable results. 

Doing so, and further assuming that the observed split- 
ting corresponds to the surface rotational frequency of 
15 CMa, we derive a rotation period of 9.35 d, which trans- 
lates into an equatorial rotational velocity of 49 ± 8kms~^ . 
The measured v sin i of 15 CMa is 40kms"^ (Uesugi & 
Fukuda 1982), making it possible that the three nonradial 
modes of the star indeed belong to the same k and £. 

For KZ Mus we confirmed the mode identifications by 
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Handler et al. (2003) . Unfortunately, the £ value of the weak- 
est known independent mode of the star can still not be 
uniquely constrained. Repeated observations in the future 
may help to solve this problem as measurements from dif- 
ferent years can bo analysed together, if homogeneous filter 
sets (particularly bluewards of the Balmer jump) are used. 
If a mode identification of /4 can be derived, we believe 
it is likely that theoretical asteroseismic investigations can 
be undertaken. Of course, additional measurements in the 
future may also help to reveal modes beyond our present 
detection limit (1.5 mmag). 

To our regret, the asteroseismic prospects of f3 CMa do 
not seem to be very high. With only three modes photomet- 
rically detected and only two identified with their £ value, 
a detailed exploration of its inner structure seems unlikely. 
Perhaps this bright star is better studied by means of spec- 
troscopy, for instance as done by Briquet et al. (2006). 

15 CMa seems more worthwhile for further study. We 
have proven the existence of a radial mode, which is of 
tremendous importance for restricting the star's location in 
the HR diagram. We have also found evidence that the re- 
maining three nonradial modes may form a rotationally split 
triplet, which can be the starting point for an asteroseismic 
study. What still needs to be done is to prove that this is 
indeed a mode triplet and to determine its £ value. 

An independent clue towards the spherical degree of 
this possible triplet may come from its frequency asymmetry 
due to the second-order effect of rotation, which is different 
for different i. Following the definition of the asymmetry 
of a frequency triplet /_ < /o < /+ by Dziembowski & 
.Jerzykiewicz (2003): 

Aobs = f-+f+ -2/(), 

we obtain A = — 0.0077±0.0015cd~^ from our measure- 
ments and A = — 0.0004±0.0011cd~^ from the archival data. 
These two values are inconsistent (perhaps due to system- 
atic frequency errors in the single-site archival data) and 
can therefore not be interpreted; more time-resolved mear- 
surements are needed. 

An unusual feature in the light curves of 15 CMa is that 
the descending branch of the strongest mode's pulsation is 
steeper than its rising branch. In fact, even a "stillstand" 
phenomenon is indicated. This has been found in only one 
P Cephei star, BW Vul, to date (e.g. see Sterken et al. 1986) 
and is usually interpreted as a shock phenomenon in the at- 
mosphere due to the vicious acceleration of the stellar ma- 
terial due to pulsation (e.g. Fokin et al. 2004). 

However, there are two differences between the pulsa- 
tions of BW Vul and 15 CMa: BW Vul pulsates in a single 
radial mode, and with much higher photometric amplitude 
than 15 CMa, whose mode showing the unusual light curve 
form is nonradial. The small photometric amplitude can be 
explained under the assumption that the strongest mode of 
15 CMa is the m — component of a rotationally split mul- 
tiplet. The possible range in the inclination of the pulsation 
axis to the line of sight may lead to heavy geometrical can- 
cellation of such a mode (e.g. see Pesnell 1985). In other 
words, this mode could have high intrinsic amplitude and 
its low photometric amplitude is just a projection effect. 

The unusual light curve form of this mode naturally 
gives rise to a harmonic frequency with high relative am- 
plitude compared to the mode causing it. Combination fre- 



quencies with abnormally high amplitudes could in some 
cases also bo caused by resonant mode coupling (see Han- 
dler et al. 2006 for a more detailed discussion), whereas other 
combination frequencies may be simple light-curve distor- 
tions caused by a nonlinear response of the flux to a sinu- 
soidal displacement due to pulsation. The sum frequencies 
we found in the amplitude spectra of KZ Mus seem to be 
such "normal" combinations. Their relative amplitudes are 
in the same range as those of the "normal" combinations of 
the /3 Cephei stars 12 Lac (Handler et al. 2006) and i/ Eri 
(Handler et al. 2004, Jerzykiewicz et al. 2005), keeping in 
mind inclinational effects. Resonant mode coupling is there- 
fore not required to explain the combination frequencies of 
KZ Mus. 

Finally, we noticed the occurrence of slow intrinsic vari- 
ability in our residual light curves for /3 and 15 CMa. Such 
variations have been reported for several other (3 Cephei 
stars (see the discussion by Handler ct al. 2006), and seem 
to be a fairly common phenomenon. Despite the existence 
of very extensive data sets, its physical cause remains un- 
known. 
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